Effect of correlated noise on quasi-lD diffusion 
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Single-file diffusion (SFD) of an infinite one-dimensional chain of interacting particles has a long- 
time mean-square displacement (MSD) oc t 1 ' 2 , independent of the type of inter-particle repulsive 
interaction. This behavior is also observed in finite-size chains, although only for certain intervals 
of time t depending on the chain length L, followed by the oc t for t — > oo, as we demonstrate 
for a closed circular chain of diffusing interacting particles. Here we show that spatial correlation 
of noise slows down SFD and can result, depending on the amount of correlated noise, in either 
subdiffusive behavior oc t a , where < a < 1/2, or even in a total suppression of diffusion (in the 
limit N — > oo). Spatial correlation can explain the subdiffusive behavior in recent SFD experiments 
in circular channels. 



PACS numbers: 05.40.-a, 66.10.C-, 83.10.Rs 

First introduced in bio-physics [l| to account for the 
transport of ions through molecular-sized channels in 
membranes, the concept of single-file diffusion (SFD), 
i.e., diffusion of particles in a narrow (quasi-lD) channel 
where mutual passage is forbidden, has_been extensively 
studied during the last few decades 0- 



Recent advances in nanotechnology has stimulated a 
growing interest in SFD, in particular, in the study of 
transport in nanopores. Thus SFD was observed in ex- 
periments on diffusion of molecules in zeolite molecular 
sieves [5j. Zeolites with unconnected parallel channels 
serve as a realization of the theoretically investigated 
one-dimensional systems. Molecular diffusion of tetraflu- 
oromethane in zeolite AIPO4-5 was studied by pulsed 
field gradient NMR spectroscopy @. The channel diam- 
eter was of the order of 0.73 nm, whereas the diameter 
of the CF4 molecules was 0.47 nm. Thus mutual pas- 
sage of molecules was excluded and it was found that the 
movement of an isolated particle, after a short ballistic 
period, was determined by the stochastic interaction with 
the channel walls yielding diffusional behavior, and the 
MSD was shown to obey an anomalous diffusion law. 

Despite the seemingly simplicity of SFD phenomenon, 
its theoretical description involves considerable difficul- 
ties caused by the necessity to solve the Fokker-Planck 
equation for the probability density distribution in a mul- 
tiparticle system. Several studies tackled this problem for 
different model systems. Thus in Ref. [7j SFD of hard- 
core particles in a one-dimensional lattice was theoreti- 
cally investigated by the method of "tagged vacancies" 
for an infinite chain, a finite chain with reflecting bound- 
aries and chains with periodic boundary conditions. It 
was shown that for an infinite chain the mean-square dis- 
placement (MSD) of a diffusing particle depends on time 
as (Ax 2 ) oc t 1 ! 2 . The same result was obtained in Ref. [|| 
and, it was demonstrated in Ref. that this result did 
not depend on the nature of the interparticle interaction. 
Note that this result was obtained assuming overdampcd 
dynamics of diffusing particles while the problem remains 



unsolved for the general case. For the intermediate case 
of neither overdamped, nor underdamped dynamics nu- 
merical simulations were presented in Ref. [9(. 

In case of a finite chain with reflecting boundaries, the 
long-time MSD was shown Q to become constant (i.e., 
independent of time t) which means a total suppression 
of diffusion. This result can be understood considering 
the fact that in a finite system (single-connected, non- 
periodic and possessing no rotational degree of freedom) , 
no matter how large, only a finite set of configurations of 
particles can be realized. Thus for a long enough observa- 
tion time, the system accesses all the states in configura- 
tion space. An exact analytical proof of this result using 
the Bethe-ansatz was recently presented in Ref. [IOJ. A 
number of numerical studies on the SFD phenomenon 
have been performed using Langevin dynamics (see, e.g ., 
Ref. [n}) including SFD on periodic substrates [l2Hl4f . 

It is hard to fulfill the SF condition and investigate 
SFD experimentally for atomic systems. Alternatively, 
experimentalists used micro-meter sized colloidal parti- 
cles in narrow channels. A well-defined SFD model was 
created fl5j by confining paramagnetic colloidal spheres 
of 3.6 /zm in a set of circular trenches (7 /xm in width and 
33 to 1608 /iin in diameter) fabricated by photolithogra- 
phy. Using video microscopy, the trajectories of individ- 
ual particles were followed over long periods of time and 
revealed a SFD behavior characterized by MSD oc i 1 / 2 . 

SFD was recently realized even on a macroscopic scale. 
The advantage of using a macroscopic system is that it al- 
lows one to easily fulfill the SF condition and observe the 
motion of particles using simple tools such as an optical 
microscope. The diffusion of such macroscopic charged 
metallic balls (of radius R = 0.4 mm and mass m = 2.15 
mg) was investigated in Ref. [16j |. The balls were inter- 
acting electrostatically and moving in a circular channel 
such that mutual passing was forbidden, while mechan- 
ical shaking induced an effective temperature. The au- 
thors found that the system of interacting balls exhibited 
a subdiffusive behavior oc i 1 / 4 , i.e., slower than that pre- 
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dieted for SFD and observed in colloidal systems. 

Here we investigate the effect of correlations of noise 
that causes diffusion. In particular, we study SFD of 
a finite-size system (i.e., interacting particles diffusing 
in a closed circular channel) driven by a mixture of a 
totally uncorrelated noise and a spatially-correlated (al- 
though uncorrelated in time) noise. Being coordinate- 
independent, spatially-correlated noise can be also re- 
ferred to as "global". We show that the presence of a 
spatially-correlated component dramatically changes the 
SFD behavior, and it can also be used as an alterna- 
tive explanation of the slowing down of diffusion observed 
in experiments with circular channels driven by artificial 
stochastic noise modeling temperature. 

The model. — Let us consider SFD of interacting par- 
ticles in a circular channel driven, by non-correlated and 
spatially- correlated stochastic forces. The fact that both 
forces are stochastic implies the absence of time correla- 
tions between the forces of the same component as well 
as of mutual correlations between different components. 
The fraction of the correlated noise to the total noise is 
given by a parameter A, where < A < 1. Therefore the 
fraction of non-correlated noise is 1 — A. The motion for 
the i-th particle is described by the Langevin equation, 



dt 2 



XF C + (1 - X)F nCji , (1) 



where r*j is the radius of z-th particle, m is its mass, 7 is 
the friction coefficient. The non-correlated part of noise 
F nc ,i = {Fn C ,i(t) , F% c i (t)} an d the spatially-correlated 
part F c = {F*(t), F^(t)} obey the following conditions, 

(F^ i (t'))=0,(F c k (t')) = 0, 

(FtA^KcA 1 ")) = WbTSuShMH - 1"), 

(F*{t')F l c {t")) =2 7 fcTfc >I 5(f -t"), 
(F*(t')Fl A (t")) = 0, (2) 

where fc, I = x,y are the coordinate indices and i, j are 
the particle indices. Such a choice of the correlation rela- 
tions ([2]) ensures that the system eventually approaches 
an equilibrium state with temperature T at any time [HJ . 
The interparticle interaction potential U(rij) in our cal- 
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culations was taken a Yukawa potential, U(r) = — - . 



Here re = 1/A where A = 4.8 • 10~ 4 m [161] . In our 
numerical calculations, other parameters of the system 
were also taken from the experiment [l6j |. The general 
conclusions are independent of the exact functional form 
of the interaction potential. The confinement poten- 
tial, U con f(ri), which restricts the radial motion of the 
particles within the circular channel is taken parabolic: 
U con f(ri) = P{r c h — ri) 2 , where (3 is chosen from the con- 
dition V gs = fir 2 ,, where ro is the radius of the particle 
and V gs is the ground state energy, i.e., the energy of 
the configuration when all the particles are distributed 
equidistant near the circular channel "bottom" . 



Long-time SFD in a finite-size chain. — In Ref. 0], 
SFD of noninteracting (i.e., hard-core) particles diffus- 
ing in a ring was studied theoretically using the "tagged 
particle" method. It was shown that the long-time MSD 
grew as 2Dt, where the diffusion coefficient D was in- 
versely proportional to the density of particles. 

We will demonstrate that this result is also justified for 
a closed (e.g., circular) chain of interacting particles. For 
this purpose, let us first notice that SFD in a finite system 
has two characteristic lengths: (i) the average interpar- 
ticle distance, a, and, as distinct from infinite systems, 
(ii) the length of the chain, L. The corresponding hierar- 
chy of times raising from the two characteristic lengths, 
implies that along with the "long-time" regime when a 
particle diffuses for a distance of the order of the average 
interparticle distance a (during time t a ), there is another 
time scale, tL, which determines the long-time behavior. 
This is the time required for a particle to diffuse over 
a distance comparable to the system length, L. Note 
that for large systems with N ^> 1 particles, t^ 3> t a 
(and tL — >• 00 when N —> 00). It is clear that due to 
the single-filing condition, diffusion of all particles will 
be approximately the same for t > tL and thus can be 
described by a "collective" coordinate R = (1/N) n. 

In case of a circular chain, one can define this collec- 
tive coordinate as an average angle of "rotation" of the 
system as a whole, <j) = ('Pi) — J2i The angu- 

lar diffusion of zth particle is described by the system of 
equations which follow from Eqs. ([!]), 



(3) 



^ = -7'^-£^xW(r y) 



dt 



+ Xn xF c + (l- X)n x F nc 



where J; is the moment of inertia of i-th particle, and 
7' = jr 2 h . Taking the average over all the particles, we 
have: 



— Ji4>i - J4> = ~i<t> + Min 
+ AM C + (1-A)M„ 



(4) 



where we assume J ~ Ji, i.e., neglect the radial disper- 
sion, and 

M mter = ~^J212 PiX VU(rij) = 0, (5) 



M r 



N, 



ddfi x F c = 0, 



channel 



M nc = £ fi X F n 



(6) 
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It is worth noting that expression ([5]) is exact due to the 
sine theorem. This well-known result from the dynam- 
ics of a conservative system is applicable to the "conser- 
vative part" of our system (although the whole system 
is not conservative). Relation ([6]) means that the field 
of stochastic spatially-correlated force is homogeneous in 
space, so the integral of this force along any closed con- 
tour is equal to zero. In other words, the action of the 
spatially-correlated force on half of particles in the chan- 
nel is on average compensated by its action on the other 
half of particles. Thus the equation of motion for the 
average angle <f> that characterizes the collective rotation 
of the system is: 



J(j> = -j4> + (1 - X)M nc . 



(8) 



This equation is isomorphic to the equation of motion of 
a free particle driven by a stochastic force. As a con- 
sequence, the long-time SFD behavior of a free particle 
can be projected on the angular diffusion described by 
Eq. (0) . Integrating Eq. ([5]) we obtain the expression for 
the long-time angular MSD: (Acj) 2 ) = ID^t. Thus we 
can conclude that long-time SFD (i.e., when t > t L ) re- 
lated to the collective rotation of the system of particles 
in a circular channel is determined only by the stochastic 
spatially-non-corrclatcd force. In other words, the in- 
terparticle interaction and stochastic spatially-correlated 
force do not influence long-time diffusion. On the other 
hand, if the system is driven only by spatially-correlated 
noise (A = I), the stochastic term in Eq. ([5]) is zero, thus 
resulting in the total suppression of long-time SFD. 

We stress that these preliminary conclusions are based 
on Eqs. ([3])-© for "collective" coordinate which are 
justified in the limit when t — > oo and N — > oo and can 
be used only for a qualitative prediction of the diffusive 
behavior in large (N ^> 1) finite systems in the long- 
time limit t >t^. Our exact results presented below are 
obtained by numerically integrating Eqs. (JT|). 

The results of numerical calculations of the MSD of 
the SFD of N particles in a circular channel are shown 
in Fig. 1 for two limiting cases: (i) for N = 20, A = 0, 
i.e., for purely non-correlated noise (Fig. 1(a)), and (ii) 
for N — 20, 30, and 45, A = 1, i.e., for fully spatially- 
correlated noise (Fig. 1(b)). In case A = 0, the initial 
fast grow (i.e., oc t 2 followed by oc t 1 ) of the MSD is fol- 
lowed by the pronounced cx t 1 / 2 "long-time" (for infinite 
systems) region (for t a < t < ti) gradually changing to 
oc t 1 asymptotic behavior for t > tt (Fig. 1(a)) which is 
common for any finite chain and any A. The short-time 
behavior for A = 1 is rather similar, while for longer time 
the t 1 growth rapidly changes to a subdiffusive behav- 
ior oc t a , with < a < 1/2, gradually approaching a 
i°-plateau for large N [T7| (see Fig. 1(b)), in agreement 
with Eq. ©. 

Spatially- correlated vs. non- correlated noise. — In our 
numerical study of angular diffusion in a circular chan- 
nel the standard deviation of an average angle (i.e., the 
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FIG. 1: (Color online) The average MSD (in log-log scale) 
of a particle diffusing in a circular channel with N particles, 
driven by: (a) purely stochastic noise (A = 0), for N — 20, and 
(b) fully spatially-correlated noise (A = 1), for N — 20, 30, 
and 45. (j> is measured in radians (also in Fig. 2). For A = 0, 
the initial fast grow oc t 2 and oc t 1 of the MSD is followed by 
the oc t 1,/2 -region (for t a < t < th) which gradually evolves 
to the long time (t > ti) oc t 1 asymptotic behavior (a). For 
A = 1 (b), the t 1 growth rapidly changes by a subdiffusive 
behavior oc t a , with < a < 1/2. 



MSD), 



N p N c 



E (< A <A - <**ij)i) ( 9 ) 



'■J 



was calculated as a function of time, t. Here (...) denotes 
averaging over all particles of the given ensemble, (•••)«, 
over various ensembles, and (...) t over time. The number 
of ensembles was chosen typically from 100 to 300, for 
N = 5 to 45 particles in the ensemble. 

Results of simulations for different A show (Fig. 2(a)) 
that when the fraction of the correlated noise increases, 
the region of fast growth (similar to the behavior of a free 
particle) of SFD oc t expands, while the region of subd- 
iffusive growth, i.e., oc t 1 ' 2 (typical for SFD of a particle 
in an infinite chain), is essentially narrowed and rapidly 
changes to even stronger subdiffusive mode oc t a , with 
< a < 1/2. Thus we observe a considerable slowing 
down of diffusion on a long-time scale when increasing 
the fraction of spatially-correlated noise. Note that for 
some value of A, i.e., A w 0.6, the long-time a becomes 
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FIG. 2: (Color online) The average MSD (in log-log scale) of 
a particle, driven by a mixture of non-correlated and spatially- 
correlated noise, for: A = 0, 0.2, 0.4, 0.5, 0.7, 0.9, and 1 (a). 
With increasing A, the region of fast growth oc t expands, 
while the region of subdiffusive growth, i.e., oc t 1 ^ 2 , shrinks 
and rapidly changes to a stronger subdiffusive mode oc t a , 
with < a < 1/2, including oc t 1/4 for A = 0.6 (b). 

Ri 0.25, i.e., leading to a pronounced region in the MSD 
oc t 1 / 4 (Fig. 2(b)). In the limit of completely spatially- 
correlated noise (A = 1), the calculated MSD (A</> 2 ) sat- 
urates for large time t (for sufficiently large number of 
particles, see Fig. 1(b)). 

The observed saturation of the MSD for long-time scale 
looks similar to that found in a finite chain where the 
motion of particles was restricted in the direction of the 
chain that in turn was related to the fact that for long 
time (of the order of the time of diffusion of a particle 
over the whole system) the trajectory of the particle filled 
practically all the configurational space accessible to the 
given particle. However, in spite of the similarity in the 
long-time SFD behavior with, e.g., a finite chain of par- 
ticles in the presence of reflecting boundaries the 
physics behind the observed slowing down of the diffu- 
sion in our system is completely different. Our system is 
periodic and it possesses a rotational degree of freedom 
which leads to a single-particle-like long-time asymptotic 



SFD behavior (as shown in Fig. 1(a)). Adding a fraction 
of spatially-correlated noise effectively suppresses the ro- 
tational degree of freedom and thus the SFD in the long- 
time limit. 

Our results provide an alternative explanation of the 
subdiffusive behavior (i.e., slower than oc t 1 / 2 ) of par- 
ticles diffusing in a circular channel driven by an ar- 
tificial stochastic force. For example, in SFD experi- 
ments with charged balls [l6j], the system temperature 
was modeled by random acoustic waves produced by a 
set of loudspeakers placed under the substrate with a cir- 
cular channel. Thus the whole sample was shaked and, 
consequently, different particles were not driven individ- 
ually. The spatial correlation of noise (although being 
uncorrelated in time) in this system arises, in our opin- 
ion, due to two factors: (i) the collective excitation of all 
the particles in the system by a limited set of sources of 
noise (note that even a large number of different sources 
would lead to spatial correlation of noise, due to the in- 
terference between the waves), and (ii) shaking the whole 
sample containing the channel with all the particles em- 
bedded in a single-file chain. Note that as was suggested 
by the authors of Ref. [l6| and confirmed by recent nu- 
merical simulations @, the interparticle interaction in 
a non-overdamped regime perhaps could also lead to a 
slowing down of the long-time SFD in circular channels. 

Similar conditions, leading to spatially-correlated 
noise, can also be realized in other artificial systems that 
exhibit single-file behavior. 

Conclusion. — We have investigated the effect of spa- 
tial correlation in noise on single-file diffusion of inter- 
acting particles diffusing in a finite-size circular chan- 
nel where the long-time diffusion is characterized by two 
time regimes, t a < t < ti and t > ti. We demon- 
strated that in the limit of purely non-correlated noise, 
the mean-square displacement characterizing the diffu- 
sion rate shows a oc t 1 / 2 (for t a < t < t L ) and a oc t 
(for t > t£) behavior independent of the type of the in- 
terparticle interaction. Adding a fraction of spatially- 
correlated noise slows down the diffusion and could re- 
sult in either subdiffusive behavior oc t a with a < 1/2 
(for t a < t < ti), or even in the total suppression of 
diffusion on a long-time scale (for N — > oc). Our model 
provides an alternative explanation of experimentally ob- 
served slowing down of the diffusion rate of interacting 
particles driven by artificial noise sources. 

We acknowledge discussions with M. Saint- Jean. This 
work was supported by the "Odysseus" program of the 
Flemish Government and FWO-V1. 
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